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Abstract We explore the possibility that the fast and exotic negative ions 
in superfluid helium are electrons bound to quantized vortex structures, the 
simplest being a ring. In the states we consider, the electron energy is only 
slightly below the conduction band minimum of bulk helium. To support our 
proposal we present two calculations. In the first, we show that the electron 
pressure on the vortex core is insufficient to cavitate the helium and form 
an electron bubble. In the second, we estimate the equilibrium radius of the 
vortex ring that would bind an electron and find it is much smaller than the 
electron bubble, about 7 A. The many exotic ions reported in experiments 
might be bound states of an electron with more complex vortex structures. 

Keywords superfluid helium • negative ions • fast ion • exotic ions 


1 Introduction 

Recent experiments by Maris and coworkers [T] have strongly substantiated 
claims of “fast” and “exotic” negative ions in superfluid helium made over 40 
years ago [211^ . Not only do these latest experiments exactly reproduce the 
relative mobilities of some 14 species first identified by Ihas and Sanders [3], 
they have resolved another four peaks in time-of-flight experiments as well 
as a signal arising from a continuously distributed mobility. That the same 
phenomena are observed with improved experimental techniques, including 
the use of carbon nanotube electrodes, argues strongly that these ions are 
intrinsic to the electron-helium system. 

The elusive nature of these negative ions, whose mobility is up to a factor 
of six greater than that of the electron bubble, is compounded by the absence 
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of a plausible model for their structure. The puzzle may be summarized as 
follows. Higher mobility implies smaller size, and suggests some entity that 
localizes a significant fraction of the electron wave function in a small volume. 
But the number of chemical species to be formed from helium atoms that 
could serve in this capacity are few and cannot account for the abundance 
of mobility peaks observed. A more speculative proposal, Maris’ fissioned 
electron bubble [S], is not a tenable model as it is based on the adiabatic 
principle of quantum mechanics in a situation where it does not apply 

This paper explores a proposal for the structures of the fast and exotic ions 
that appears to have been overlooked. Consider how these ions are created 
in the laboratory. Electrons are injected into liquid helium after they have 
diffused through the vapor above the liquid surface. Based on the energy 
parameter of the Meyer-Reif distribution [5] of the diffusing electron in a 
field, Wei et al. estimate the electron energy in the range 3-6 eV at the time 
of injection [T]. In the liquid the electron creates excitations, slowly losing 
energy and increasing its wavelength. When the energy has dropped to a 
particular value, Vb ~ 1 eV, the electron stops diffusing because below this 
energy the electron has no propagating states. This cutoff energy is usually 
referred to as the conduction band minimum of liquid helium. 

In the adiabatic approximation, where the electron motion is so fast that 
helium atoms may be treated as static (in a configuration typical of the 
liquid), the constant Vb corresponds to the mobility edge of Anderson local¬ 
ization. This interpretation of Vb is better suited to our discussion, because 
the nature of the localized state accessed by the electron, even while the 
helium atoms are effectively static, determines the eventual structure of the 
ion. 

The simplest geometry of the localized electron state is a roughly spherical 
wave function nestled in one of the larger cavities of the fluid. The normal 
electron bubble is formed, a self-trapped state, when on longer time scales the 
walls of this cavity are expanded by the zero-point pressure of the spherical 
electron wave function. Our proposal for the fast and exotic ions starts with 
the observation that the development of the localized state might be more 
complex and have alternate end-products. 

Because the electron is much lighter than a helium atom, energy loss is 
a slow process and passage through the mobility edge is not abrupt. Hydro- 
dynamical evolution of the fluid therefore does not switch on at a particular 
time, but develops by degree as the character of the electron wave function 
becomes increasingly localized. The hydrodynamical forcing by the electron 
can impart impulse to the fluid, a condition the superfluid can support only 
by being threaded with vorticity. Considering the fact that the energy of a 
10 A radius vortex ring is only about 0.01 eV, the production of vorticity can 
be significant even for electrons that have lost just a fraction of the mobility 
edge energy of Vb. 

The helium density at a vortex core is diminished and provides a line of 
weak attraction to the electron. And since in two dimensions an arbitrarily 
weak attractive potential binds a quantum particle, a straight vortex will 
confine the electron in its transverse plane. It follows that a localized distri¬ 
bution of vorticity of sufficient size, say a ring, will form a bound state with 
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the electron. The electron energy for these states is only slightly below Vq 
because there is no large cavity that holds most of the electron wave function. 
We speculate that the fast and exotic negative ions may be explained by this 
class of bound state. 

When modeling the final equilibrium state of the electron in the presence 
of vortices we encounter the following conundrum. In the Anderson localiza¬ 
tion picture, the electron sees the helium atoms as quenched disorder and 
has decaying wave function 

tZ/oc (1) 

where the correlation length 

c oc (Fo - E)-^ (2) 

has exponent v Ki 1 This behavior, of the envelope of the wave function, 
applies to all the different realizations of “disordered” helium atoms that 
arise in the fluid. On the other hand, in the conduction band picture used 
extensively in the past, bulk helium is treated as a continuum medium that 
imposes a uniform potential Vq on the wave function. The decay of the wave 
function will then have exponent v = 1/2. We have not tried to resolve this 
discrepancy and will adhere to the potential function model. 

The best available modeling for the electron-helium interaction is based 
on density functionals [?]■ This paper uses more elementary modeling and 
serves to motivate such a study. We carry out two calculations. In the first 
we show that an electron weakly bound to a linear vortex will not cavitate 
the core to transform it to a normal electron bubble. In the second, we find 
the equilibrium size of a vortex ring that binds an electron. Although the 
crudeness of our modeling of the vortex core introduces uncertainties, the 
resulting size is small and makes the ring a compelling candidate for the fast 
ion. 

More complex vortex structures - links and knots - look like they might 
explain the multiplicity of ions observed in experiment. No calculations have 
been made in this paper to support that hypothesis. A class of structures 
that holds promise m is the different ways a ring with multiple quanta of 
circulation may resolve itself into multiple filaments, possibly knotted and 
linked, all singly quantized. Owing to their parallel-oriented vorticity, such 
vortex tangles are less likely to reconnect and might support stable periodic 
motion. If the exotic ions corresponded to dynamic equilibrium structures 
such as this, they would be examples of true perpetual motion. 

We have not explored extensions of the vortex-bound electron states that 
could explain the recently observed [1] continuous component of the mobil¬ 
ity distribution. However, the fissioning wave function proposal [I] suffers 
from the same adiabaticity-breakdown concerns as its original application 
to photo-excited electron bubbles [S] . In brief, Maris and coworkers consider 
the scenario where an electron entering the helium with energy exceeding Vq 
becomes a superposition of reflected and transmitted waves, with amplitudes 
determined as in the textbook problem of a quantum particle encountering 
a fixed potential step. The transmitted part of the wave function then cav¬ 
itates the helium to produce a bubble of variable size, consistent with the 
transmission amplitude. This analysis fails to acknowledge the fact that, in 
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general, the electron-helium system must be described by a joint wave func¬ 
tion. When the electron coordinate in this joint wave function is above the 
helium surface, the helium atoms are not compelled to preserve the bubble 
geometry that prevails when the electron coordinate is under the surface. 
The correct superposition is thus always a combination of normal electron 
bubble, and, a reflected-wave electron and cavity-free helium. 


2 Calculational framework and definition of scales 


We model the helium as a continuum fluid of uniform mass density p that 
forms a sharp interface at the empty cores of vortices. The interfacial surface 
energy, of a vortex ring of radius R and core radius a, is 

Hs = An'^Raa, (3) 


where a is the surface tension and we have assumed R ^ a. The energy 
of the helium in this model is just the kinetic energy in the corresponding 
classical flow field for vortex circulation h/M {M is the helium mass): 


= ^ (log(8i?/a) - 2). 


(4) 


The electron Hamiltonian is 


It o 

He = - 

2m 


H(r) 


(5) 


where m is the electron mass and the potential V is zero at points within 
distance a of the vortex core and Vq otherwise. We will let Vq be our unit 
of energy and define our unit of length oq in terms of the associated kinetic 
energy: 


K) = 




277100 


( 6 ) 


For Vb = 1 eV, the scale oq = 1.95 A is comparable to the radius of a vortex 
core. Scaling the surface and flow energies in terms of our energy and length 
scales, we obtain the following two dimensionless parameters: 

= 0.00326 (7) 


a = 477 ^ 


Vb 




= 0.000878. 


( 8 ) 


With our choice of energy and length scales the total Hamiltonian of the 
electron-vortex-ring system is 


H = He+He+Hh (9) 

= -V^ -f u(r) -k aRa + PR (log(8i?/a) - 2), (10) 

where u(r) = 1 except inside the tube of radius a surrounding the ring, where 
it is zero. Although the magnitudes of the surface and flow terms are small 
due to their overall scales, we will see in the next Section that the ground 
state energy of He is comparably small for the case of interest. 
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3 Stability of the vortex core 

In cylindrical coordinates appropriate to the ring geometry, the electron 
Hamiltonian has the form 


He = Ho + Hi, 


( 11 ) 


where 


f d'^ 9^ \ 

^o = -(^ + ^j+-(D.) 

(12) 

H ^ ^ 

^ rdr' 

(13) 


In this Section our interest is the core stability of a large ring where we can 
neglect the term Hi which accounts for the curvature of the ring. This term 
is restored in the next Section where we determine the equilibrium radius of 
the ring. 

In terms of the distance 


s = ^{r- i ?)2 + z 2 


(14) 


from the (effectively straight) vortex core, the ground state electron wave 
function is 


'f'(s) oc 


Jo{ks), s < a 
Ko{ks), s > a, 


(15) 


with eigenvalue 


Hq = = 1 — k^. 


(16) 


The value of k is determined by matching the logarithmic derivative at s = a. 
Defining the variables x = ka and y = na, the system of equations to be 
solved is 


Jijx) _ Ki{y) 
^Jo(x) ^ Koiy) 

= a^, 


(17) 

(18) 


whose solutions are functions of the single parameter a. The electron energy 
is given by 


Ho = 



(19) 


In the limit of small a, that is, cores smaller than the scale oq ~ 2 A, 
equation Cll) can be approximated as 


2 log (f 2 /) ’ 


( 20 ) 


where 7 is Euler’s 
mation we find 


constant. Solving the system of equations in this approxi- 



1 - (4/a2)e-4/“'-2^, 


( 21 ) 









6 




a (A) a (A) 


Fig. 1 Excess energy above the conduction band minimum of the electron-vortex- 
ring bound state as a function of core radius for two sizes of ring: = 5 A (left) 
and = 20 A (right). 


showing that the energy of the bound state (below the conduction band 
minimum) is exponentially small in the quantity 1/a^. 

In Figure 1 we show the behavior of the total energy H — Vq as a function 
of the core radius a for two values of R. The energy has a local maximum 
for a of order 1 A, and decreases weakly for small a and strongly for large 
a. For large a the vortex core is expanded by the electron pressure, the fluid 
undergoes cavitation, and an electron bubble is formed. We are interested 
in the case of small a, where the electron pressure is exponentially feeble 
and the equilibrium structure is determined by the energetics of just the 
helium (balance of surface tension and Bernoulli pressures). If we ignore the 
electronic contribution entirely, the equilibrium core radius is given by 

a = /3/a = 0.269, (22) 

or O. 53 A in physical units. An electron that finds itself in the vicinity of a 
vortex that is already small could not use it as a cavitation nucleus. The 
thermal activation barrier, for cavitation by uniform expansion of the core, 
is about 10 K for a 5 A ring and about five times this value for a 20 A 
ring. Although it has negligible effect on rings with small cores, an electron 
would nevertheless be bound when the ring is large enough. The effect of 
the bound electron on the size of the vortex ring is the task we take up in 
the next Section. In these calculations we treat the core radius as a fixed 
phenomenological parameter. 


4 Electron-vortex-ring equilibrinm 

In the previous Section we saw that binding of an electron to a vortex lowers 
its energy, below the conduction band minimum, by an amount that behaves 
exponentially with the core radius. This was in the limit of a large ring, where 
the transverse extent of the electron wave function is small compared with 
the ring’s radius of curvature. Since the electron wave function falls off expo¬ 
nentially, the effects of ring curvature on the binding energy is exponential 
as well. Finding the lower limit on the ring size, below which it no longer can 
bind an electron, is therefore a delicate mathematical exercise. 
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Two corrections must be applied to the calculation of Section [S] The first 
is to recognize that for the quasi-2D Hamiltonian Hq, an image vortex-core 
must be included at r = — i? in the r-z plane in order to satisfy the boundary 
condition 


d'l' 

dr 


= 0 . 

r—0 


(23) 


Writing the center of the true vortex (in the r-z plane) as Ri and the image 
as R 2 , the electron wave function outside the core is 


tf'out(r) =A{Ko{K\r - Ri|) -h Ko{K\r - R 2 I)), (24) 


where r is a point in the r-z plane and A is the normalization constant. 
Since most of the wave function is outside the core in the case of interest, we 
determine A by ignoring the contribution from the interior wave function. It 
will prove to be very convenient to express all the calculations in this Section 
in terms of the dimensionless parameter 

q = kR. (25) 

The normalization condition takes the form 


1 = A^R^fiq), 

(26) 

where 

pOO poo 

X 2 

f(q)=4TrJ dy J xdx {^Ko{qy/ix - y'^)-\-Koiqy/{x-\-1)'^y'^)j . 

(27) 

The second correction in the binding energy calculation is the term Hi of 
the electron Hamiltonian ((1311. Treating this term as a first-order perturbation 
to be applied to the unperturbed wave function , we obtain the correction 

.00 .00 / IQ^\ 

(28) 

= Trj dz'R^^^ir = 0,z) 

(29) 

g{q) 1 

/(g) 

(30) 

where 

pCO 

g{q) = 87r dy -\-y"^). 

Jo 

(31) 


Combining the first two orders in the computation of the electron energy, 
we obtain 

g(g) \ 1 


H, = Ho + (Hi) = l+[-q 


f{q)J i?2(g)- 


(32) 


We will see shortly how the ring radius R is also a function of q, owing 
to wave function matching at the core radius. The sign of the multiplier of 
1/R^ determines whether the electron is bound. Numerically we find there 










is binding for q > 0.38. However, the true bound on q is smaller because we 
know the second order correction of a ground state energy (not calculated 
here) is always negative. 

The wave function matching condition, in the case of a small core, is the 
same as equation (EOl) but with an extra term arising from the image core: 


(fca)^ —1 

2 “ \og{^Ka) - Ko{2KRy 


(33) 


Using kR = q and na = q{a/R), and solving for R we obtain 


^g2/(fea) 

2 


(34) 


The electron wave vector in the core, k, depends weakly on R through the 
relation 


= 1 — = \ 


from which we get the expansion 

g2/(fea)= ^ g2/a2 



- id/Rr, 

(35) 

{Ray ^ ) 

(36) 


For the range of ring sizes and the parameter q we will need, the correction 
terms are negligible and we can use the approximation 


R{q) 


y 




x/5 


qe 


-Ko{2q) 


(37) 

(38) 


where we have introduced a small parameter associated with the electron 
energy: 



(39) 


The electron binding energy is now expressed in terms of this scale and a 
function of q: 




o2ifo(29) 


(40) 


The possibility of the electron stabilizing a vortex ring against shrinking 
down to zero radius depends critically on the parameter 5. By (1391) the latter 
depends sensitively on the dimensionless core radius a, and to a lesser extent 
on the conduction band minimum Vq through the definition of the length 
scale oq ([S]). In the most optimistic scenario, when both of these are large, 
<5 is largest and is best able to compete with the flow and surface tension 
energies of the vortex. For example, if we take a = 1.5 A and Vq = 1.3 eV, 
then S = 0.009. This number is reduced by about one order of magnitude if 
the core radius is decreased to 1.2 A. 

Figure 2 shows the electron binding energy (go]) vs. ring radius (l38ll plot¬ 
ted parametrically as functions of q and scaled by their respective powers of 
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Fig. 2 Scaled electron binding energy (HOI) and corresponding scaled vortex ring 
radius dSll) plotted as functions of the parameter q. The electron is bound to the 
ring for q > 0.38 and R > 0.208/\/5. 


5. Although the details are more subtle than the case of the normal bubble, 
the electron zero-point energy favors a large ring. It is this weak lowering 
of the energy with ring size that competes against the helium energy, which 
we take to be the energy of the flow in a ring vortex o, where the core¬ 
surface energy contribution is included by adjusting the core radius a. The 
flow energy depends weakly on a as well as Vb, through the scale oq and 
the parameter /3 We choose a = 1 A and Vb = 1 eV, giving /3 Ri 0.0009. 
Figure 3 shows the total electron-vortex-ring energy (relative to Vb) vs. ring 
radius for S = 0.006. The equilibrium value of the radius is R = 6.7 A. In 
Table I we give binding energies, radii, and the parameter q for equilibrium 
rings at other values of 6. These results show that when an equilibrium bound 
state exists, the vortex ring has radius 6.7 A, and that probably no such state 
exists when <5 is below 0.002, because then q is below the value for binding. 

The electron energy and the hydrodynamic energy of the flowing helium 
are sensitive to very different characteristics of the vortex core. To model the 
former we introduced S, and treated it as an independent parameter even 
though it is nominally related to the hollow-core radius a of the hydrodynamic 
model through equation (1551) . To make the best case for the electron-vortex 
bound state we selected optimal values for these parameters. Better modeling 
by density functional methods will eliminate this freedom and provide a more 
critical test of the proposed ion model. 
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Fig. 3 Total energy of the electron-vortex-ring bound state as a function of ring 
radius for the case S — 0.006. Equilibrium radii for various S are given in Table 1. 


5 Creation and transport of electron-vortex-ring bound states 

There can be no doubt that an electron injected into liquid helium above 
the conduction band minimum can be captured as a bound state to a large 
enough vortex ring (or something topologically more complex). As long as 
the vortex core radius is below some minimum, of the same order as a free 
vortex, the electron pressure is insufficient to cavitate the core and form an 
electron bubble. What happens after capture is less certain. 

Because the total energy of the electron-vortex-ring bound state (calcu¬ 
lated in Section 11 is an increasing function of radius for large radius, over 
time the ring will decrease its radius through emission and scattering of 
phonons and rotons. The most intriguing scenario is that this process will 
cease, with equilibrium established when the ring has radius of order 7 A. 
This possibility rests on the reliability of our bound state energy estimate, 
which unfortunately scales with the small parameter 6 whose value is sensi- 


Table 1 Properties of equilibrium electron-vortex-ring bound states for various 
values of the parameter 5. 


5 

H-Vo (K) 

R{k) 

<7 

0.010 

-20.2 

6.72 

0.515 

0.008 

8.3 

6.74 

0.482 

0.006 

36.8 

6.73 

0.440 

0.004 

65.0 

6.66 

0.388 

0.002 

90.8 

6.38 

0.310 
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lively dependent on the core radius and the conduction band minimum, Vq. 
When 6 falls below about 0.002, the electron becomes unbound {He > Vq) 
and the ring is free to decay. 

The scenario of no equilibrium electron-vortex-ring state has trouble ex¬ 
plaining the transport measurements of the fast and exotic ions. There are 
two cases to be considered. For large electric field E, when the electric force 
is balanced by a drag force proportional both to the radius R and velocity v 
of the ring, we get the relations 


E oc Rv, 


(41) 


u oclog i? oc Fie (42) 

R 

where we have used the velocity u of a ring advected by its own flow field. In 
the weak field case, when the drag must be modeled as stochastic, the ring 
can decay before its impulse and size are restored by the electric force. Now 
the electron is liberated and free to be recaptured by vorticity elsewhere in 
the helium. Thus, while a decaying-exponential drift velocity is consistent 
with the high field measurements of Eden and McClintock m, the low field 
mobility would appear to lack the robust characteristics recently established 
by Maris and coworkers [T]. 

The most promising scenario then, of electron-vortex-bound states as the 
elusive fast and exotic ions, is for the electron to establish an equilibrium 
structure. These would be metastable states, as the electron energy is nearly 
1 eV above the energy of the ground state electron bubble. Nevertheless, their 
stability is made possible by two things. First, because the electron is in the 
ground state with respect to the potential function of the vortex-threaded 
helium, the state cannot decay radiatively. Second, decay via tunneling to an 
expanded core that can cavitate is suppressed by the relatively large size of 
the equilibrium ring. 

The starting point for discussing the transport properties of an equilib¬ 
rium ring is to realize that such a ring would have zero velocity (in the 
absence of a field) even while its impulse p is non-zero. By general principles. 


dH _ dH dR 
dp dR dp 


(43) 


since the total energy is stationary at the equilibrium ring radius. The cartoon 
of flow fields in Figure 4 shows how this can be understood in terms of force 
balance. The two terms in H contribute opposing generalized forces: 


dH dHh , dHe 

'hr ~ 


(44) 


The first force is the inward Magnus force per unit length of vortex times the 
length of the ring. This is balanced by the outward force produced by the 
bound electron. 

Like the ordinary electron bubble, the zero-velocity electron-ring bound 
state acquires a finite drift velocity in the presence of a field. Two features 
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Fig. 4 Contrasting flow fields around a free vortex ring (left), and a ring that 
binds an electron (right). The free ring is advected by its own flow field and moves 
to the right with velocity v as shown. Becanse the core is at rest with respect to 
the flow, the flow aronnd it is nniform and there are no Bernoulli forces. The flow 
pattern around the u = 0 ring that binds an electron, on the right, is obtained 
by transforming the free ring on the left to the frame moving with velocity v. 
This results in a flow asymmetry around the core and unequal Bernoulli pressures 
that produce a net inward Magnus force F. Equilibrium of the core is maintained 
because the bound electron generates a net outward force of equal magnitude. 


suggest the drag on the ring, due to phonon-roton scattering, would be sig¬ 
nificantly smaller for the ring than the electron bubble: (i) its size (radius) 
is smaller, (ii) the “helium contrast” is much smaller (empty core vs. empty 
bubble). On the other hand, unlike the electron bubble, the vortex ring is 
surrounded by a non-zero flow field even in the limit of zero velocity. This 
has the effect of enhancing the effective scattering cross section. 


6 Conclusion 

It is known that vortices can bind electron bubbles. The bound states con¬ 
sidered in this paper are qualitatively different in that the electron leaves a 
much smaller footprint in the helium. As a possible model of the fast negative 
ion, we studied the binding of an electron to a vortex ring, where the effect 
of the electron is simply to stabilize the ring against decay when its radius 
is about 7 A. The helium density in this proposed ion is much more uniform 
than when the electron forms a bubble. Although the vortex ring’s flow field 
extends the range over which it scatters phonons and rotons, the small size 
and more uniform density of the proposed ion would increase its mobility 
significantly over that of the electron bubble. 
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The most perplexing experimental fact that any model of the exotic ions 
must address is their great number. Although we have analyzed in detail only 
the vortex ring, the same stabilization mechanism should apply to topologi¬ 
cally more complex vortex structures, such as knots and links. As in the case 
of the ring, the electron’s energy would be only slightly below the conduction 
band minimum. Even so, the energetics of the electron state is such that it fa¬ 
vors straight over curved vortex filaments, the net effect being to expand the 
vortex “tangle”. The density functional methodology [7], suitably adapted, 
might be able to test the validity of this mechanism and if confirmed, identify 
the most stable structures. 
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